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Abstract 



A microscopic theory of zero wavelength (q = 0) interaction in finite-size 

systems is proposed. Its exact solution interpolates between the Coulomb 

blockade and the perturbative Altshuler-Aronov theory, in the strong and 

weak interaction limits respectively. The tunneling density of states and the 

quasiparticle life-time are calculated. The physical nature of the (q = 0) 

component of the interaction is discussed. 
PACS numbers: 72.10.Bg, 73.50. 
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Much of the activity in the field of mesoscopic physics has focused on zero-dimensional 
(d = 0) weakly disordered systems, where the relevant energy scales turn out to be smaller 
than the Thouless energy, E c = HD/L 2 , D being the diffusivity and L the system's linear size 
(for a review see Ref. Studies of transport as well as of thermodynamic properties reveal 
the importance of electron-electron interactions || . The role of the latter and the interplay 
with disorder have been underlined by the perturbation theory developed by Altshuler and 
Aronov (AA) who consequently explained the zero-bias anomaly in tunneling into 

a macroscopic conductor. Some ten years ago the observability of single electron effects in 
small capacitance devices has been proposed |5]||, leading to strong suppression of tunneling 
at low external bias. Nazarov [0] has pointed out to the common physics between the 
perturbative AA anomaly and the non-perturbative Coulomb blockade. 

There are several interesting issues pertaining to d = systems. One is the tunneling 
density of states (DOS). To study this quantity it is possible to employ perturbation theory 
in the spirit of AA; in doing so one should recall that the relevant wavenumbers are quan- 
tized with a special role played by the q = component of the interaction. This component, 
pertaining to fluctuations of the global charge, is unscreened. When the interaction associ- 
ated with this component is strong the phenomenological Coulomb blockade theory may be 
employed, but it is also desirable to have a microscopic theory at hand. Another interesting 
issue involves the electron inelastic life-time. 

Here we present an exactly solvable microscopic model for interacting electrons in zero 
dimensional systems. Our model interpolates between the Coulomb blockade limit and 
the zero-bias anomaly. The results for the latter turn out to be in line with a simple 
generalization of the AA calculations [|J. We thus find the tunneling DOS as function of 
energy, temperature and coupling strength. We also evaluate the single particle life-time, 
associated with the width of the single particle levels. Finally we discuss the nature of the 
q = interaction term, V(q = 0). Relating it to the system's capacitance we are able to 
establish a connection with previous studies of current fluctuations through an ultrasmall 
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capacitance junction |8|||. 

We consider a system described by the Hamiltonian 



(1) 



where {e a } represent exact single particle eigenenergies; the parameter N depends on the 
positive background and the external fields (e.g. gate voltages), and needs not to be an 



integer. The imaginary time single particle Green function is given by [TU 
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where Z(fi) is the partition function, given by the same integral without pre-exponential 
factors. Next we perform Hubbard-Stratonovich transformation by means of the imaginary 
Bose field u(t), yielding 
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with the same transformations in Z(fi). Here Z^(/j,) and Q^(ti,t 2 , /i) are the partition 
function and the Green function, respectively, in the presence of a time dependent spa- 
tially uniform potential, cr(r), superimposed on the single particle problem |lTJ. Introduc- 
ing the Matsubara representation (for the boson field uj m = 2iiimT, m = 0,±1,±2...), 
c"m = Jo dra(r) exp{cj m r} , we obtain 

zW(fi) = ZM(v-<To) 

(n, r 2 , = gm (ri , T2 , „ - a )Jn d ^-°°\ (4) 

where Z^{jj) and (?M(e n ,/i) are partition function and Green function of non-interacting 
electrons. Substituting in Eq.([3]), we have 



-^ m -(exp{— u m 72}- exp{-w m Ti})] 
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Here fiM is the (non-interacting) thermodynamic potential. The integral over the static 
component, o"o, describes the smooth transition between the grand canonical ensemble with 
the chemical potential fi (V = 0) and the canonical ensemble with a given number of 
particles [No] (V = oo). For sufficiently large systems (or high temperature), ^ + y ^> (3 
(A -1 = — d 2 Vt^ /djj? - is the mean level spacing), differences between the two ensembles are 
negligible and the integral over <7o may be evaluated within a saddle point approximation (at 
o"o = (To). The m/0 components are readily evaluated. The single particle Green function 
assumes the form 

Gain - r 2 , fi) = G [ a ] (n - r 2 , /x - W )T>(ji - r 2 ), (6) 

where the Green function of the auxiliary field a, representing the Coulomb boson, is given 
by 

V{t) = e~^( W "^); \t\<(3 (7) 

(it is periodic with a period 0). In the Matsubara representation the Green function is 
superimposed with the Coulomb boson, see Fig. [I], to yield 

G a {e n ,n) = T^2G [ a ] (e n - u m ,-p)V(uJ m ), (8) 
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where /! = // — Wo, e n = 2iri(n + 1/2)T and D(w m ) is evaluated as the Fourier transform of 
X>(t) (in imaginary time), over the period [0,0]. The related spectral function is defined by 
B(uj) = 2QV a (uj) [|T2J (T> A is the advanced Green function of the Coulomb boson) and is 
given by 



/ 27T / _{u+V/2f „(^-v/2) 2 \ 

B{u) = sl— [e 2vt -e 2vt y (9) 

Employing the spectral representation of G&, A%\ the total Green function assumes the 
following Lehmann representation 

_ . . Iff 00 de' du', #\ a rni / / N coth^ + tanh^ 



We are now in a position to evaluate the DOS, v(e) = —n 1 Z)a^^a( e )- Analytically 
continuing Q a , we finally find 

where is the DOS related to Q$. For non-interacting electrons (with or without disor- 
der), i/[°l(e) ~ const, substituting Eq. (§) into Eq. ([[]]), we obtain 
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Note that this result is independent of the disorder in the system. We have thus derived an 
expression for the DOS which interpolates between the perturbative zero bias anomaly and 
the semiclassical Coulomb blockade. In the limit of weak interaction (V <C T) this reduces 
to 

= i - ^cosh" 2 ^. (13) 
z/M AT 2T y ' 

This turns out to be the zero-bias anomaly of AA if one substitutes d = in their analysis 
H. For strong interaction (V ^> T) 



'2ttT v , e 



Ke) _ I \/ir e ~" cosh 2T ; e « ^« ^> (14) 
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The exponential suppression of the tunneling DOS is a direct manifestation of the Coulomb 
blockade. In the extreme case of zero temperature one has 



v(e) T -^^6(\e\--). (15) 

To discuss the quasiparticle lifetime due to global charge fluctuations, we need an explicit 
form of the electron Green function. To this end we substitute the spectral function of non- 
interacting electrons, A$(e, ~p) = 2ir5(e — e a +~p), into Eq. fllPp , perform integrations and 
do analytical continuation (e n — > e + i8). The result is 
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where /(e) is the Fermi function, x = e — e a + ~p + i5 and w(z) = exp{— z 2 }erfc(— iz) is 
the error function [P^j . We note that the w(z) does not have poles at finite energies. We 
nevertheless may consider the width of the peak of the imaginary part of Q R as a measure 
of the quasiparticle life-time. As a result 
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Both the DOS and r qp are measurable in tunneling experiment []. 

Finally we discuss the nature of V = V(q = 0). Consider the equivalent circuit shown 
in Fig. ||. Its equilibrium noise spectrum (T = 0) is given by the fluctuation-dissipation 



(UtatU tot ) = iu (-L + Z(lu) ) . (18) 



theorem 



\iujC 

The voltage fluctuations on the system (capacitor) are Uq = U tot /(l + iuCZ), so the corre- 
sponding correlation function is 

%)^(W^ U ^ H . (19) 

We may now add to our original Hamiltonian, Eq. ([[]), a Gaussian noise term elicit) Y^ a a tt a a 
(the correlator of the noise is given by Eq. fll9|)). Integrating over the environment noise 
and following through the above derivation of Q a , we shall have to replace V — > V — K(u m ) 
in Eq. (|^). It is natural to identify V = e 2 /C, where C is the total capacitance of the 
d = system. In this case our noise-modified Coulomb boson (after Gaussian integration 
over a m in Eq. @) reproduces the result of Refs. and for fluctuations in ultrasmall 
capacitance junctions. 

We are indebted to A. Altland and A. M. Finkelstein for the useful suggestions. This 
research was supported by the German-Israel Foundation (GIF) and the U.S.-Israel Bina- 
tional Science Foundation (BSF). 



1 Note that the AA result for the dephasing time Q scales as ^ oc T 2 /( 4 ~ d ) D^ d /^~ d \ For d = 
it resembles our Eq. (|l7|). 
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Note added after this work has been completed : in a recent preprint by Levitov and 
Shytov (SISSA N 9501130), these authors present an effective semiclassical action through 
which they study tunneling into a d = 2 system. Their q = limit agrees with our exact 
result. 
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FIGURES 

FIG. 1. The single particle Green function in the presence of the Coulomb boson, T>. 

FIG. 2. An equivalent circuit, consisting of the system's capacitance C and a noisy environ- 
ment. 
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